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Abstract. In the interstellar medium the turbulence is believed to be forced mostly through 
supernova explosions. In a first approximation these flows can be written as a gradient of a 
potential being thus devoid of vorticity. There are several mechanisms that could lead to vorticity 
generation, like viscosity and baroclinic terms, rotation, shear and magnetic fields, but it is not 
clear how effective they are, neither is it clear whether the vorticity is essential in determining 
the turbulent diffusion acting in the ISM. Here we present a study of the role of rotation, shear 
and baroclinicity in the generation of vorticity in the ISM. 
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The study of the interstellar medium (ISM) is strictly connected with that of super- 
i— - n ' novae explosions, being this one of the most important phenomena in determining its 
, dynamic. Such events act on scales up to ~ 100 pc injecting enough energy to sustain 

■ turbulence with velocities of ~ 10 km/s. An understanding of the generation of turbu- 
lence and vorticity in the ISM is a basic requirement in order to formulate a theory 

■ for the production of interstellar magnetic fields. The first step to perform numerical 
simulations of these explosion is to assume the presence of pure potential forcing. In- 

O ' deed a supernova explosion can be seen at first glance as a pure spherical expansion 
£j that is driving turbulence in the ISM. In our simulations we use the Pencil Code, 
| http :// pencil- co de . googl ecod e . com/I We simulate spherical expansions following the work 

i i of Mcc & Brandenburg] (|2006l) . We analyze flows that are only weakly supersonic and use 

a constant and uniform viscosity in an unstratified medium. We solve the Navier-Stokes 
equations in the presence of viscosity and with a potential forcing V</> where 4> is given 
; by 

^" '• (j){x,t) =0o Nexp{[x-Xf(t)] 2 /R 2 } . (0.1) 

Here x = (x,y,z) is the position vector, Xf(t) is the random forcing position, R is the 
_* radius of the Gaussian, and A is a normalization factor. We consider two forms for 
the time dependence of Xf. First, we take Xf such that the forcing is (^-correlated in 
time. Second, we include a forcing time 5tf mco that defines the interval during which Xf 
t-H . remains constant. We study then the effect of rotation on such environment. Rotation 
causes the action of the Coriolis force, 21~i x u, in the evolution equation for the velocity. 
We investigate flows with Reynolds numbers up to 150. We find that there is a clear 
production of small scale vorticity when we use both a finite and 5-correlated forcing. 
In this last case we observe some spurious vorticity also at very low Coriolis number, 
due probably to some numerical artifact. We then consider the effects of shear on the 
potential flow. In the presence of linear shear with u s = (0, Sx, 0), the evolution equation 
for the departure from the mean shear attains additional terms, — it s V • u — u • Vu s . 
Also in this case we observe the production of vorticity but we do not find clear results 
for small values of the shear. We have thus observed a production of vorticity due to 
rotation and shear, but these effects are not able to pro duce enough vort icity under 
the physical condition of the ISM, like those described by iBeck et al.l (|l996h . When we 
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Figure 1. Images of temperature, entropy, baroclinic term (VT x Vs) z and normalized vertical 
vorticity for a two-dimensional run with Stf OICC c s o/R = 0.1 at an instant shortly before the 
second expansion wave is launched (top row), and shortly after the second expansion wave is 
launched (second and third row). In the second and third row the vorticity production from the 
baroclinic term is clear, while in the top row, ( VT x Vs) z and io z are just at the noise level of the 
calculation. Shock surfaces are well localized and the zones of maxi mum production of vorticity 
are th ose in which the fronts encounter each other. Adapted from iDel Sordo fc Brandenburg 
il20l"0T) . 

relax the isothermal condition the system is then also under the action of a baroclinic 
term. This is proportional to the cross product of the gradients of pressure and density 
and emerges when taking the curl of the pressure gradient term, p^ 1 Vp. The results we 
obtain are shown in Figure [TJ 

It turns then out that the baroclinic term is more efficient in the production of the vor- 
ticity. Moreover the biggest amount of vorticity is observed when shock fronts encounter 
each other. 

Spe aking about possible dynamo action, as pointed out by iBrandenburg fc Del Sordo 
(2009), the presence of vorticity does not seem to affect the diffusion of magnetic fields 
differently than a complete irrotational turbulence. Nevertheless the vorticity plays an 
important role in dynamo processes so it is still important to address the problem of the 
generation of the vorticity investigating the possible role of other effects. 
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